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We present the first official release of the nCTEQ nuclear parton distribution functions with errors. 

The main addition to the previous nCTEQ PDEs is the introduction of PDE uncertainties based 
on the Hessian method. Another important addition is the inclusion of pion production data from 
RHIC that give us a handle on constraining the gluon PDE This contribution summarizes our 
results from [1] and concentrates on the comparison with other groups providing nuclear parton 
distributions. 
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1. Introduction 

Nucleons and nuclei can be described using the language of parton distribution functions 
(PDFs) which is based on factorization theorems [2-5]. The case of a free proton is extremely 
well studied. Several global analyses of free proton PDFs, based on an ever growing set of pre¬ 
cise experimental data and on next-to-next-to-leading order theoretical predictions, are regularly 
updated and maintained [6-10]. The structure of a nucleus can be effectively parametrized in 
terms of protons bound inside a nucleus and described by nuclear PDFs (nPDFs). These nPDFs 
contain effects on proton structure coming from the strong interactions between the nucleons in a 
nucleus. Similarly to the PDFs of free protons, nuclear PDFs are obtained by fitting experimental 
data including deep inelastic scattering on nuclei and nuclear collision experiments. Moreover, 
as the nuclear effects are clearly dependent on the number of nucleons, experimental data from 
scattering on multiple nuclei must be considered. In contrast to the free proton PDFs where quark 
distributions for most flavors together with the gluon distribution are reliably determined over a 
large kinematic range, nuclear PDFs precision is not comparable due to the lack of accuracy of the 
current relevant data. In addition, the non-trivial dependence of nuclear effects on the number of 
nucleons requires a large data set involving several different nuclei. Nevertheless, nuclear PDFs 
are a crucial ingredient in predictions for high energy collisions involving nuclear targets, such as 
the lead collisions performed at the LHC. 

In this contribution we present the new nCTEQ nuclear PDFs that were recently released [I] 
and compare them with analyses from other groups providing nPDFs [11-13]. All the details of 
the analysis can be found in ref. [1] here we will mostly concentrate on the differences with other 
nPDFs. 

2. nCTEQ global analysis 

In the presented nCTEQ analysis we use mostly charged lepton deep inelastic scattering (DIS) 
and Drell-Yan process (DY) data that provide respectively 616 and 92 data points. Additionally we 
include pion production data from RHIC (32 data points) that have potential to constrain the gluon 
PDF. To better asses the impact of the pion data on our analysis two fits are discussed: (i) the main 
nCTEQlS fit using all aforementioned data, and (ii) nCTEQ15-np fit which does not include the 
pion data. The framework of the current analysis, including parameterization, fitting procedure and 
precise prescription for the Hessian method used to estimate PDF uncertainties is defined in ref. [1] 
and we refer reader to this paper for details. 

In both presented fits we use 16 free parameters to describe the nPDFs, that comprise 7 gluon, 4 
M-valence, 3 <i-valence and 2 <7-|-« parameters. In addition, in the nCTEQl 5 case the normalization 
of the pion data sets is fitted which adds two more free parameters. Both our fits, nCTEQlS and 
nCTEQ15-np describe the data very well. Indeed, the quality of the fits as measured by the values 
of the ;t;^/dof (0.85 and 0.87 for the nCTEQlS and nCTEQ15-np fits respectively), confirms it. 
Figure 1 shows the bound proton PDFs resulting from the two fits. The nuclear correction factors 
(left panel) clearly show that the pion data impact the gluon distribution, and to a lesser extent the 
My, dv and s PDFs. Also, the inclusion of the pion data decreases the lead gluon PDF at larger 
X (> 10^^), and increases it at smaller x whereas the error bands are reduced in the intermediate 
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Figure 1: Comparison of the nCTEQlS fit (blue) with the nCTEQ15-np fit without pion data 
(gray). On the left we show nuclear modification factors defined as ratios of proton PDFs bound 
in lead to the corresponding free proton PDFs, and on the right we show the actual bound proton 
PDFs for lead. In both cases scale is equal to 2 = 10 GeV. 


to larger x range. For most of the other PDF flavors, the change in the central value is minimal 
(except for a few cases at high-x where the magnitude of the PDF is small). For these other PDFs, 
the inclusion of the pion data generally decreases the size of the error band. 

The description of the fitted data by the nCTEQl 5 fit can be seen in Fig. 2 where we display 
ratios of F 2 structure functions for different nuclei as well as the corresponding data. Finally, in 
Fig.3, we display the comparison for the ratios of: F 2 structure functions for iron over deuteron 
(left), and pion yields in DAu over pp collisions at RHIC (right). In both cases the nCTEQlS fit 
provides a very good description of the experimental data. 

3. Comparison with other nPDFs 

We now compare the nCTEQlS PDFs with other recent nuclear parton distributions in the 
literature, in particular HKN07 [11], EPS09 [12], and DSSZ [13]. Our data set selection and 
technical aspects of our analysis are closer to that of EPS09 on which we focus our comparison 
in the following. Our results for the nuclear modifications of the lead PDEs as well as the bound 
proton lead PDEs themselves are shown in Eig. 4 for different flavors at the scale 2 = 10 GeV. 
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Figure 2: Comparison of the nCTEQl 5 NLO theory predictions for R = /F^' with nuclear target 

data. The bands show the uncertainty from the nuclear PDFs. 


For most flavors, d, s and g, there is a reasonable agreement between predictions from 
different groups. In particular, for the gluon, there is a larger spread in the predictions form the 
various PDF sets; we can see a distinct shape predicted by the nCTEQlS and EPS09 fits whereas 
HKN07 and DSSZ have similar, much flatter behavior in the small to intermediate .r region and 
deviates form each other in the higher .r region; however, all these differences are nearly contained 
within the PDF uncertainty bands. 

On the other hand, examining the u- and cf-valence distributions one can see a very different 
pattern. Three of the PDF sets, HKN07, EPS09 and DSSZ, are very similar across the x range, 
while the nCTEQlS set displays a marked deviation from this behavior.* At small-x (< 10^*) 
the Uy is above the others while the dy is below; and, for large-x the exact opposite behavior is 
observed. This trend persists across all Q values. In essence, the average value of the Uy and dy 
nuclear corrections are comparable to the other groups, but individually the corrections are very 
different. 

This difference highlights an essential feature of the nCTEQlS fit; namely, that the Uy and dy 
are allowed to be independent, whereas other groups assume the corresponding nuclear corrections 
to be identical. Certainly there is no physical motivation to assume the Uy and dy nuclear corrections 
to be universal; we exploit this freedom and fit the two PDEs independently. To demonstrate that 
the difference in the shape of the valence distributions is not an artifact of our methodology, but 

'This difference is significantly reduced when we consider full nuclear PDFs, + (A — Z)/Zf"l^, 

which are the quantities entering cross section calculations (see Fig. 25 in [1]). 
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Figure 3: (a) Ratio of the F 2 structure functions for iron and deuteron calculated with the nCTEQlS 
fit at = 5 GeV^ overlaid with fitted data [14-18] and results from EPS09 and HKN07. (b) 
Comparison of the nCTEQlS and EPS09 fits with the PHENIX [19] and STAR [20] data for the 

’ Fa^^^/dpTcly 


ratio = 


rather the result of the data and the physically motivated freedom in valence nuclear corrections, 
we performed an additional fit, nCTEQ15-mod, in which the Uy and dy parameters have been 
tied together in order to obtained a universal correction. This result is shown in Eig. 5 where the 
nCTEQ15-mod fit aligns closely with the EPS09 prediction. The of the nCTEQ15-mod fit 
increases to 677 (;^^/dof = 0.94) as compared to 611 (;^^/dof = 0.85) for the nCTEQlS fit. This 
is more than our tolerance criteria, = 35, used for defining fhe error PDEs, and indicafes fhaf 
fhe dafa prefer fhis addifional freedom. 

booking af Eig. 3 if is also clear fhaf fhe differences in fhe valence disfribufions do nol influence 
fhe abilify fo describe fhe currenfly available dafa. We can see fhaf observables in Eig. 3 are very 
well described by bofh nCTEQlS and EPS09 fifs despile fhe very differenl valence PDEs.^ This 
facl highlighfs fhe need for more dafa fo consfrain fhe nuclear PDEs. 

4. Summary 

We have presenfed fhe nCTEQl 5 sef of nuclear PDEs including a sef of error nPDEs generated 
using fhe Hessian mefhod. In addifion fo fhe sfandard DIS and DY data sets, pion production data 
from RHIC were also included, providing an additional handle on the gluon PDF 

We have compared the results of the nCTEQ15 fit with nPDFs from HKN07, EPS09, and 
DSSZ groups. While there are similarities between the nCTEQ15 fit and the other sets at a macro 
level, there are significant differences in the details. In particular, the most notable and important 
difference is the treatment of the m- valence and r/-valence nuclear corrections. While other groups 

^The very notable difference of nCTEQl 5 and EPS09 error band in Fig. 3b is caused by the fact that EPS09 include 
pion data with weight of 20 (whereas we use weight of 1) which effectively shrinks the errors. 
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Figure 4: Comparison of the nCTEQlS fit (blue) with results from other groups: EPS09 [12] 
(green), DSSZ [13] (orange), HKN07 [11] (red). The left panel shows nuclear modification factors 
for lead, and the right panel the actual PDFs of a proton bound in lead. The scale is 2 = 10 GeV. 



Figure 5: Bound proton and dy PDFs for lead for a modified nCTEQ15-mod fif using universal 
nuclear correction for both valence distributions (purple). For comparison we show distributions 
for the nCTEQlS fit (blue) and EPS09 fit (green). The scale is 2 = 10 GeV. 


use a universal nuclear correction for the valence distributions we believe that there is no physical 
reason for doing so and we treat them as independent. This additional freedom leads to an improved 
fit with significantly lower which in turn is preferred by the available data. However, to clearly 
answer this question, more data allowing the separation between u and d quarks will be needed. 
One example of such a data is the W jZ boson production in pPb collisions at the EHC, e.g. [21]. 
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Our PDFs are publicly available at the nCTEQ website: www. ncteq. org. We provide grids 
in the internal CTEQ PDS format (together with the corresponding interface) and also in the new 
LHAPDF6 format [22]. 
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